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In this work, we describe a method to realize 3D Weyl semimetal by coupling multilayers of 
honeycomb optical lattice in the presence of a pair of Raman lasers. The Raman lasers render each 
isolated honeycomb layer a Chern insulator. With finite interlayer coupling, the bulk gap of the 
system closes at certain out-of-plane momenta due to Raman assisted tunnelling and result in the 
Weyl semimetal phase. Using experimentally relevant parameters, we show that both one and two 
pairs of Weyl points can be realized by tuning the interlayer coupling strength. We suggest that 
Landau-Zener tunnelling can be used to detect Weyl points and show that the transition probability 
increases dramatically when Weyl point emerges. The realization of chiral anomaly by using a 
magnetic field gradient is also discussed. 


I. INTRODUCTION 

The search for topological phases that are fully gapped 
in the bulk and possess topologically protected gapless 
surface states has been an important subject in physics. 
Several new topological phases, such as 2D and 3D topo¬ 
logical insulators mm and quantum Anomalous Hall in¬ 
sulators [3], have been realized experimentally in con¬ 
densed matter systems. The gapless surface states in 
these systems are protected by certain symmetries and 
the bulk gap mm- Nevertheless, the surface states can 
be removed once the bulk gap in the system is closed. 

Interestingly, it is well known that Weyl semimetals nr 
[8], which contain nodal points, the Weyl points, in the en¬ 
ergy spectrum also support topologically protected gap¬ 
less Fermi arcs on the surface. A pair of Weyl points in 
a Weyl semimetal can be regarded as a pair of topologi¬ 
cal defects in 3D momentum space with opposite chirality 
and possess linear energy spectrum near the Weyl points. 
The Fermi arcs, which connect the Weyl points in the sur¬ 
face Brillouin zone, cannot be removed unless the Weyl 
points are merged together by strong perturbations mm- 
Due to the presence of the Weyl points in the bulk and the 
gapless Fermi arcs on the surface, Weyl semimetals are 
expected to exhibit interesting novel phenomena which 
are absent in fully gapped topological phases such as the 
chiral anomaly IS HO] and unconventional quantum os¬ 
cillations m- Recently the experimental realization of 
Weyl semimetal phase is under widely studying in real 
materials [HIIS]. 

On the other hand, great advancement has been made 
in realizing topological phases with cold atoms in optical 
lattices. Using Raman laser coupling pAll2Q] and opti¬ 
cal shaking lattice [2T1123] . Chern insulators have been 
successfully realized experimentally [24l |25] using both 
schemes. Methods on realizing 3D topological insulators 
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in optical lattices have also been proposed [MIEI]. In 
this work, we outline a practical scheme to realize and 
detect Weyl semimetal phase in optical lattices. 

We base our proposal on the already realized honey¬ 
comb optical lattice [28l |29] and show that by introduc¬ 
ing non-trivial hopping between the neighboring honey¬ 
comb lattice layers, it is possible to generate either one 
or two pairs of Weyl points in the system. What is im¬ 
portant in our condtion is that the momentum transfer 
Ap from Raman beams possesses nonzero components 
both along and perpendicular to the honeycomb lattice 
planes. Within each individual lattice layer, in-plane mo¬ 
mentum transfer generates non-trivial Peierls phases and 
can drive each layer into a 2D Chern insulator phase. Im¬ 
portantly, when individual layers are coupled with non¬ 
trivial phases generated by the Raman beams, the Chern 
insulators can become trivial for certain out-of-plane mo¬ 
mentum kz. Since the bulk gap has to be closed for 
topological phase transitions from the topologically non¬ 
trivial Chern insulators to trivial insulators as a function 
of out-of-plane momenta kz, there are Weyl points in the 
momentum space separating the trivial and non-trivial 
regimes. The number of Weyl points, corresponding to 
the number of the band-gap closing points as a function 
of kz, is tunable by changing, for example, the Raman 
detuning and the interlayer coupling. We suggest that 
Landau-Zener transitions can be used to detect the emer¬ 
gence of the Weyl points. We also discuss how Landau 
levels and the chiral anomaly can be realized by intro¬ 
ducing a magnetic field gradient in the 2 :-direction. 


II. THE SCHEME 

We first consider a multilayer honeycomb optical lat¬ 
tice generated by laser beams [28l [29] . 

V = —Vx cos^{qx 0/2) — Vx cos^(gx) — Vy cos^ (qy) 

- 2a^/VxVY cos{qx) cos{qy) - cos^{qJ), (1) 
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where Vx , Vx are intensities of lasers X, X and 

Y, pointing along x, —x and ^-directions, respectively, q 
is the wave vector of the lasers. When Vx ^ Vy ^ hz, a 
honeycomb lattice is created with two lattice sites A and 
5 in a unit cell, as shown in Figj^. The energy offset A 
between A and B sublattices can be tuned through the 
phase difference 0 between X and X lasers. The stand¬ 
ing wave Z, with wave vector Qz, stacks the honeycomb 
optical lattice as seen in Fig{^, and controls the inter¬ 
layer coupling strength. The energy offset A is tuned to 
be large enough such that the direct tunnelling between 
neighboring A and B sites is greatly suppressed in the 
absence of Raman beams. 

A pair of Raman beams with frequency detuning 6uj = 
uji — co ’2 ~ A/fi are added to restore the resonant tun¬ 
nelling between neighboring intralayer sites, as illustrated 
in Fig{^. The momentum transfer Ap = ki — k 2 as¬ 
sociated with the resonant tunnelling introduces Peierls 
phases exp [iAp • (r + r')/2] between intralayer neigh¬ 
boring sites r and r'. This introduces synthetic mag¬ 
netic flux into the optical lattices nnsni as depicted in 
Figigi. Using Floquet analysis and a unitary transfor¬ 
mation which eliminates the spatial dependence of the 
Peierls phases (for details, see Appendix [A|), the effective 
Hamiltonian can be written as 

+ h.c.. (2) 

Here, i, j labels the position of lattice sites within a layer 
and I is the layer index, annihilates an atom 

at sublattice A (B) of site i and layer 1. which satis¬ 
fies due to the Raman fields |3T], denotes 

the intralayer nearest neighbour hopping amplitude and 
is the next nearest neighbour hopping amplitude 
between A (B) sublattice sites, and are the inter¬ 
layer hopping amplitudes and e = ^ {hduj — A), originat¬ 
ing from the Raman detuning, can be regarded as the 
on-site energy difference between the sublattices. The 
Peierls phases (^ij and are associated with intralayer 
and interlayer hopping, generated by the in-plane and the 
out-of-plane components of Ap, respectively. Explicitly, 
one finds 

4>ij = ■ (ri,i - i-i,;), = ^Ap • (rj,i - rj,i+i). (3) 

In particular, (j)ij is responsible for making each isolated 
layer of honeycomb lattices a Chern insulator and (j)x^ 
which originates from the phase difference between ad¬ 
jacent layers due to nonzero vertical momentum trans¬ 
ferred in the Raman process, induces unconventional in¬ 
terlayer coupling which is the key for generating the Weyl 
semimetal phase. The accumulation of Peierls phases 
along closed paths gives rise to synthetic magnetic flux 
patterns as shown in Fig{^ and e. 



FIG. 1: (color online) (a) AA-stack of multilayer honeycomb 
optical lattice, (b) The intralayer optical potential. The big¬ 
ger dark blue and the smaller light blue pockets represents 
the A and B sites respectively, (c) The Raman laser assisted 
resonant tunneling between nearest neighbors. The recoil mo¬ 
mentum carried by photons generates Peierls phases in the 
hopping amplitudes, (d) The synthetic magnetic flux of 
each triangle resulting from the Peierls phases in the a:y-plane. 
(e) The synthetic magnetic flux of each rectangle in the yz- 
plane. 


III. THE WEYL POINTS 

To illustrate the emergence of Weyl points, we first 
consider an ideal version of experimental honeycomb lat¬ 
tice (brickwall lattice). For simplicity, we set 1 = t and 
(in the following NNN hopping ti is along vi or 
V 2 while t 2 is along vi — V 2 ). We use the coordinate sys¬ 
tem established as shown in Figj^ for convenience and 
the momentum transfer is set to be Ap = {PxiPyiPz)- 
If we further label the A{B) sublattices as correspond¬ 
ing to a pseudo-spin cr, then the Hamiltonian takes the 
following form in momentum space 

H (k) = do (k) / + d (k) • o- (4) 

with 

do(k) = 2ti cos kxa cos -pxCt + 2ti cos kya cos -Pya 

+ 2^2 cos{kx + ky)a cos ^ {px + Py)a 
+ 2t±_ cos kzd cos (j )±_, 

d^(k) = tcosi {kx - ky)a, 

dy(k) = —2t sin - {kx + ky) a + t sin -^{^x ~ ky)a^ 

dz(}Y) = 2ti sin kxa sin -pxa + 2ti sin kya sin -pya 

+ 2t2 sm{kx + ky)asm^{px Y Py)a 

— 2t±_ sin kzd sin (j)±_ -h e. 


( 5 ) 
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FIG. 2: (color online) (a) The Chern number C of an isolated 
layer with t 2 — 3ti = 0.3t as a function of the Zeeman term m 
and recoil momentum Ap = (0,^2/, 0). The purple dots denote 
m{kz = 0) = e, while the dashed lines correspond to the range 
of m for general kz when interlayer coupling is turned on. (b) 
The energy difference between the upper and the lower bands, 
using parameters from a realistic optical lattice, is shown. It 
is evident that two pairs of Weyl points locate in the plane of 
ky = kx — TV. (c) Energy spectrum in the ky = kx — tt plane 
with e = 0 where two pairs of Weyl points emerge, (d) One 
pair of Weyl points emerges with e = 0.2t. In (c) and (d) we 
take ti = O.lt, t2 = 0.3t, t± = 0.26t and Ap = (0,7r/a, 7r/d) 
for calculation. 


where a and d are lattice constants in the x^-plane and 
along i-axis. A dependent Zeeman term m = e — 
2t^^mkzd^m(pj_ appears in the effective Zeeman fields 
along ^-direction. For fixed m, the Chern number 

C = ^ y* dkxdkydiy) • {dxdCk) x dxd{k.)) (6) 

can be defined as the topological invariant of the system 
with fixed kz and (i(k) = d(k)/|d(k)|. 

The phase diagram of an effective 2D system with 
as a function of fixed m and Py is sum¬ 
marized in Figj^, where the colored region has non-zero 
Chern number. When the interlayer coupling t± is zero, 
m = e with pya G (0, 27r) U (27r, 47r) makes \dz\ >0 for all 
kz. As a result, the system is fully gapped in the bulk. 

As the interlayer coupling is increased to a finite value, 
m as well as the bulk gap change as a function of kz. 
From Eqj^ it is clear that the Weyl points appear when 
{kxa, kyo) = (—l^r, —|7r) or kyo) = (|7r, —^tt) and 
at a value of kz where dz{kz) = 0 such that |d| = 0. For 
systems with specific t± and Ap, the ranges of m{kz) 
as a function of kz, where kz ranges from —ir/d to Tr/d, 
are schematically indicated as vertical dashed line (i) for 
e = 0 and (ii) for e = 0.2t in Figjit. The number of 
Weyl points is determined by the number of crossings of 



FIG. 3: (color online) (a) The two planes which contain the 
Weyl points are depicted. The color represents the energy 
difference between upper band and lower band and indicate 
the position of the Weyl points. Driven in the (110) direction, 
atoms at the bottom of the lower band will move along the 
arrows in the k space and penetrate the two planes. Atoms 
in a unit Bloch oscillation encounter two local band gap min¬ 
ima and undergo two Landau Zener tunneling events, (b) 
Three typical cases of Landau-Zener transition. The path in 
the upper and bottom panel contains Weyl points with differ¬ 
ent locations, while that in the middle panel is fully gapped. 
The dashed arrow denotes that the Landau-Zener transition 
probability is small due to the existence of the energy gap. 


the dashed lines and the topological phase boundary in 
Figj^. The crossing is controllable by tuning the inter¬ 
layer coupling, the detuning and the momentum transfer 
Ap in the Raman process, so the Weyl semimetal phase 
with either one or two pairs of Weyl points becomes avail¬ 
able. The energy spectrum for systems which contain one 
or two pairs of Weyl points are depicted in Figj^ and d, 
respectively. The range of the Zeeman terms m as a func¬ 
tion of kz in Figj2fc and d correspond to the vertical lines 
(i) and (ii) in Figj^ respectively. 

In the following, we consider a realistic optical poten¬ 
tial with 0 = and V^x,x,y] ~ [7,0.5,2]F^i^, where 
Eji = 12m is the recoil energy [28l [29]. The hop¬ 
ping parameters of the optical lattices are calculated 
from the maximally localized Wannier function method 
(see Fig ]A3| in the Appendix). In realistic band struc¬ 
tures with the Raman laser strength Vq = Eji, two 
pairs of Weyl points emerge successively when the am¬ 
plitude of is increased. Using t± = 0.09Eji^ Ap = 
(0, 27r/a, Tr/d), e = 0, we plot the energy gap between the 
lower and the upper bands in Figj^. Here the energy 
gap is plotted for the planes kz = 0, kz = tt, ky = kx — x 
and ky = —kxE kxi + kyi with i = 1,2. The four crossing 
points {kxi,kyi,0) and {kxi,kyi,7r) between the orthog¬ 
onal planes are the Weyl points. It is evident that the 
four band touching points represent the Weyl points. It 
is important to note that, using the realistic parameters, 
there is a wide range of experimentally accessible param¬ 
eter space of Eji and Vz in which the Weyl semimetal 
phase can be realized. 
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IV. DETECTION OF WEYL SEMIMETAL 
PHASE 

Our proposed scheme for realizing Weyl semimetal con¬ 
sists of two steps: 1) create isolated 2D Chern insulators; 
and 2) couple the isolated layers with non-trivial hop¬ 
pings. To confirm the emergence of the Weyl phase, we 
shall first detect the Chern insulator and then ascertain 
the closing of bulk gap when interlayer tunneling is in¬ 
creased. As we show below, both can be carried out with 
Bloch oscillations. 

It is important to note that the nonzero Chern num¬ 
ber comes from the integral of Berry curvature over the 
first Brillouin zone. Thus one way to reveal the topolog¬ 
ical order is to measure the Berry curvature, as was pro¬ 
posed theoretically in [32j |33] , and performed experimen¬ 
tally [ISl |29] . In the Bloch oscillation induced by a con¬ 
stant force, atom cloud, in addition to the group velocity 
parallel to the force, would obtain an anomalous veloc¬ 
ity that is transverse to the force due to nonzero Berry 
curvature m- After one full Bloch cycle, this Berry cur¬ 
vature generated anomalous velocity would slightly drive 
the atom cloud transversely. This transverse drift, which 
is measurable experimentally, can manifest the topologi¬ 
cal order of Chern insulators. 

After the confirmation of 2D Chern insulator phase for 
the isolated layers, the interlayer coupling can be gradu¬ 
ally enhanced through modulating the intensity of laser 
Z to close the bulk band gap [29] . Such band gap closing 
can be monitored by the fraction of atoms in the excited 
band after a cycle of Bloch oscillations, as was done for 
isolated 2D honeycomb lattice [28l |35] . Below, we extend 
this technique to Weyl semimetal case and calculate the 
Landau-Zener probability. 

More specifically, we use the model described by H (k) 
in Eqj^ and consider an external linear potential ap¬ 
plied along (110) direction where there are two pairs of 
Weyl points, see Figj^. The total probability for inter¬ 
band transition will come from two independent Landau 
Zener events as schematically depicted in FigJ^. Since 
the transition probability in each Landau-Zener event in¬ 
creases exponentially as the gap decreases, the emergence 
of Weyl points is accompanied by a dramatic increase in 
transferred fraction when bulk band gap closes. For non¬ 
interacting fermionic atoms, as was performed in recent 
experiment [28|, the Landau Zener probability must be 
averaged over the initial distribution of atoms [35] . 

Let us denote the energy of the lowest band as F(k) 
and the energy difference between the lowest and the first 
excited bands close to band touching point as ^(k). Then 
within local density approximation, the Landau-Zener 
transition probability can be written as 

/B(k,r)<Bp^’(k)^"kd3r 

^ r d^kd^r ’ ^ ^ 

JE{\i,r)<EF ^ ^ 

where Ep is the Fermi energy and F(k, r) = F(k) + V(r) 
with F(r) the external harmonic confining potential. 



FIG. 4: (color online) (a) Landau-Zener transition probabili¬ 
ties as a function of t± for different Fermi energies. A transi¬ 
tion peak signals the emergence of Weyl points at the critical 
value of t±. The peak separates the fully gaped Chern insula¬ 
tor phase from the Weyl semimetal phase, (b) Landau-Zener 
transition probability as a function of both t± and p with 
Fermi energy Ep = 0.05t. In calculations, we take F = 0.01, 
ti = O.lt, t 2 = 0.3t, e = 0 and Ap = (0,7r/a,7r/d). 


p(k) is the Landau-Zener probability for k-state, 
7r(5f(k)-|/^ r 




AvjE 


1 — exp 


4vjF 


), ( 8 ) 


where labels the corresponding Landau-Zener events 
along (110) direction and uqj) is the velocity close to 
Weyl points. E is the strength of the linear potential. 

The transition probability as a function of for sys¬ 
tems with different Ep^ is shown in Fig|^. The momen¬ 
tum transfer due to the Raman field are Pya = it and 
Pzd = IT such that the system can be a Weyl semimetal 
when t±_ is larger than some critical values. Initially, 
atoms with energy lower than Ep are occupied near the 
k = (0,0, ±7r/2(i). For small the system is fully 
gapped and the Landau-Zener transition probability is 
exponentially small. At certain critical t_L, Weyl points 
emerge on the line along the (110) direction. Due to 
the gap closing, the inter-band transition probability in¬ 
creases dramatically as shown in Figj^. However, by 
further increasing the Weyl points move away from 
\kzd\ = 7r/2 plane. As a result, the atoms driven in the 
(110) direction will not be able to access the Weyl points 
and the transition amplitude decreases. Therefore, there 
is a peak in the transition probability as a function of t± 
and the peak location indicates the critical value of t± at 
which the Weyl points emerge. In other words, the tran¬ 
sition peak separates the fully gapped Chern insulator 
phase from the nodal Weyl semimetal phase. 

The density plot of the transition probabilities as a 
function of t±_ and Py for Fermi energy Ep = 0.05t is 
shown in Figj^. The Landau-Zener transition peak lo¬ 
cates around tj_ = ^ |sin|p^a(ti — ^ 2 )], which is de¬ 
noted as the dashed line in FigJ^. As in Figj^, the 
peak separates the Chern insulator phase from the Weyl 
semimetal phase. Similar results can be obtained when 
realistic parameters are used. 
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FIG. 5: (color online) (a) Schematic illustration of uniform 
magnetic flux $ generated by a magnetic field gradient in the 
^-direction and the Raman laser assisted tunnelings. This 
flux $ coexist with which is responsible for the emergence 
of the Weyl points, (b) The energy spectrum of the Weyl 
semimetal in the presence of The in-gap gapless modes 
(denoted by the red lines) connect the gapped Weyl points 
with opposite chirality. The in-gap gapless modes denotes 
the energy spectrum formed by the zeroth Landau levels. 


V. LANDAU LEVELS AND CHIRAL ANOMALY 


direction, as it is schematic shown in Figj^. 

In the presence of this uniform magnetic field along the 
x-direction, energy spectrum formed by Landau levels are 
generated. The zeroth Landau levels result in gapless 
chiral modes along the direct ion at the Weyl points. 
The chiralities of the gapless chiral modes are determined 
by the chiralities of the Weyl points m- The energy 
levels in the presence of the synthetic magnetic flux as a 
function of are depicted in Figj^. The in-gap modes 
represent the energy of the zeroth Landau levels. 

With a constant force applied parallel to the magnetic 
field along the x-direction, atoms can be adiabatically 
pumped from one Weyl point to another Weyl point 
with opposite chirality. This is the manifestation of chi¬ 
ral anomaly. The chiral anomaly induced atom number 
imbalance between different Weyl points is expected to 
cause asymmetric quasi-momentum distribution in the 
time-of-flight measurements j37H39] . Furthermore the re¬ 
cent advancement in cold atom quantum transport mea¬ 
surement j4Ql442] can enable the study of chiral anomaly 
induced negative magnetoresistance m in atomic sys¬ 
tems. 


It is important to note that the Weyl semimetal 
phase discussed above can be used to realize the chi¬ 
ral anomaly. Considering a low-energy Weyl Hamilto¬ 
nian without any velocity anisotropy for simplicity, in 
the presence of an external magnetic field B = (H^,, 0, 0), 
we can choose the Landau gauge A = (0,0, eBy) to ob¬ 
tain iLw = X'^F [^kx<Jx — + {Skz + eBy) az], where 

X = ±1 represents the chirality of the Weyl point. 
Along the direction parallel to the magnetic field, the 
infinite Land au bands form gapped spectrum En = 
hvFSgn{n) y/2 |n| eB/h + 5k‘^, with n = ±1,±2,... and 
also gapless chiral modes Eq = X'er^kx. The gapless 
chiral modes induced by the zeroth Landau levels have 
opposite group velocity for the states at Weyl points of 
opposite chirality. When a constant force (gradient scalar 
field) is applied parallel to the magnetic field, states near 
one Weyl point can be driven to the other Weyl point of 
opposite chirality through the gapless chiral modes. In 
this case the particle number is not conserved in a single 
Weyl point but the total particle number over the two 
Weyl points remains a constant. In the background of 
the magnetic field and gradient scalar field, this viola¬ 
tion of particle number in a single Weyl point of specific 
chirality is refered as chiral anomaly [36]. In the follow¬ 
ing, we will show that a uniform synthetic magnetic field 
can be introduced in our Weyl semimetal phase to form 
the chiral Landau bands and further simulate the chiral 
anomaly. 

In our AA stacked honeycomb optical lattice, when 
a linear magnetic field gradient is induced along the z- 
direction and together with the laser assisted tunneling, 
an extra phase ^ = Ap • Vij^i = E j(l)y + is added 
to the 2 : direction hopping [19]. In our choice of 
Raman lasers, (0^,, (py^(pz) = (O 5 f 5 'tt), the extra phase ^ 
results in a uniform synthetic magnetic field along the x 


VI. CONCLUSION 

In conclusion, we introduced an experimentally feasi¬ 
ble way for realising Weyl semimetals by applying Ra¬ 
man lasers to coupled multilayer honeycomb optical lat¬ 
tices. We suggest that the measurements of Landau- 
Zener tunneling probabilities can be used to detect the 
Weyl semimetal phase and the possible exploration of chi¬ 
ral anomaly is discussed. Our scheme will lay the founda¬ 
tion for the realisation of even more exotic nodal topolog¬ 
ical phases such as Weyl superconductors [H] [45] when 
attractive interactions between atoms are introduced. 
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Note Added: At the finishing stage of this work, we 
noted that a scheme to realize Weyl semimetal phase by 
stacking 2D Harper systems in cubic lattices was pro¬ 
posed [46] . 


Appendix A: TIGHT BINDING MODEL 

Coupled by Raman field, the periodic drive system in¬ 
volving the resonant modulation can be described by a 
time-independent effective Hamiltonian [30] [31]. With 
hopping terms indicated in FiglAll the Hamiltonian 
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FIG. Al: The hopping terms considered in the honeycomb 
lattice (brickwall lattice). 


reads 


m,n,l 


+ 

+ 


2Vo sin (Ap • - 6ujt) - - 

2Vb sin (Ap • - 6ujt) + ^ 




^m,n,rrn,n,l 


+ ^02 

+ toiCilri,n,l^'>^,n,l + +h.C.^ (Al) 


where i = mvi + nv 2 + Id and = mvi + V 2 + 

Id + — Rq. Considering the gauge transformation 


U = exp f -iSujt 


m,n,l 

+ exp ( - - 


.t 




(-Ii(5wt^ ^L,n, 


f^bm,n,k 5 


(A2) 


we have 


H' = U^HU -ihU^^ 
dt 

= Ho + 


(A3) 


with 




m,n,l 


+ t 2 Kn-l,n+l,firn,n,l + t±bl^^^i_^bm,n,l 

m,n,l 

+ t02 (^bl^_^^^ j^am,n,k + 

AT/oe K^P'^rr.,n,z 2 )a^,n,z 

+ (A5) 


and 


^-1 ~ 

m,n,l 

+ to2 + A,n-l,fe^m>«>fe) 

+ (A6) 

Here we consider the effective Hamiltonian to the first 


order 


ft.=-f/o+im.-f/-j+o(y, 


then 


- u 

^eff 


' 0 , 


(A7) 


(AS) 


^ ^ (^01 + 2 ^ 02 ) ^]n,n,Z^7n,n,Z 
m,n,l 


/ /1 5 71 5 6 

A ^ 01^02 

+ ^02^]n,n,Z^^+l»^-l^^ “ (^01 + ^^ 02 ) ^ln,n,l^^:n,l 

toito 2 A 


-^2 t 

^02^m,n,Z^7n+l,n-l,Z 


nL,IL,L I ^ 5 " 

_|_ P/noT/^^’^ P * 2 ^P’(^m,n,Z Tr^_ ]^ ^ ^) t b a i 
Hr ZlQ2yrn-l,n^ ^ ^ ^ CLrn,n,rrn-l,n,l 

TPtnol/"'’^ p^4^P-(C,n,Z+rm-l,n,z)/)t 

^zro 2 r^_l^^e ’ ^ o^_^^^ iam,n,l 
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FIG. A2: The Fermi arc on the (100) surfaces of Weyl 
semimetals, (a) using the parameters corresponding to Fig.2c 
in the main text, (b) using the parameters corresponding to 
Fig.2d in the main text. 


where = Vo sin ^Ap • (r^ „ ^ - r^,-It is clear 

that the Raman lasers will introduce the Peierls phase 
exp [iAp • (r + r') / 2 ]. For the ideal model in the main- 
text, we assume = 1 2 = ^ 2 , and ig¬ 

nore the difference in AA and BB sublattice hopping 
in the effective Hamiltonian. For Ap = (O^Py^Pz), 

I sin Ap - -r^-,z)l has the same value for the nearest 

neighbor hopping, so the total Hamiltonian is simplified 
as 


H = 






+ ^ h.c.. (AlO) 


The spatial dependence of the Peierls phases is eliminated 
through the unitary transformation as follows 


Appendix B: WEYL POINTS IN REALISTIC 
OPTICAL POTENTIAL 


ai,i = ai^i exp j 

h,i = exp 

After the gauge transformation, the new Hamiltonian 
takes the form 

Vi),I 

+eE +E 

-i-h.c.. (B 2 ) 



FIG. A3: The band structure and its corresponding maxi¬ 
mally localized Wannier function for the lattice potential with 
kx,x,v = [7, 0.5,2] Er and 0 = |7r. The hopping terms ob¬ 
tained from maximally localized Wannier function makes the 
energy spectrum agree well with that from plane wave results, 
k is along kx — ky. ^ — V^cl is the wavelength of the 

in-plane laser. 


Further labeling A {B) sublattices as corresponding to a 
pseudo-spin a, then the Hamiltonian is transformed in k 
space as follows 

H (k) = = do (k) J + d (k) . < 7 . (B3) 

Weyl points will emerge when the conditions dx (k) = 
dy (k) = dz{k) = 0 are satisfied. For the case |t^ | = 
t, dx (k) = dy (k) = 0 occurs when {kxa^ kyo) = 
(—Itt, — Itt) or {kxa^kyo) = (| 7 r, — ^tt). Thus the crit¬ 
ical condition for the emergence of Weyl points can be 
determined by judging the existence of real roots for 
(i ;2 (—Itt, — Itt,= 0 or (—^tt, | 7 r,= 0 with 
given e and 0 _l.P hysically the existence of Weyl points 
corresponds to a transition from topological nontrivial 
state to topological trivial state when atoms move along 
kz direction. The Fermi arc evolution is shown in Fig |A 2 | 
Realistic optical lattice potential with Vjxxv] ~ 

[7, 0.5, 2] and ^ = | 7 r is considered to calculate the 
band structure for isolated layer and its correspond¬ 
ing Maximally localized Wannier function [47]. The 
Wannier function W (x, y, z) = w (x, y) w (z) is con¬ 
trolled by in-plane and out of plane opti cal potential 
respectively, and w {x,y) is shown in Fig |A3| with in- 
tralayer hopping parameters [toi,to 2 ,^ 2 ] = 
10-2 [-8.29,-1.65,0,0.45,1.86,-0.21,3.22] and the 
onsite potential A = —1.2Er. In the resonant case 
A = Scj, we can select the wavelength for the z direction 
standing wave to have Pzd = ir and make Vb = Er. Such 
optical potential can generate Weyl points as is shown in 
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